NACA TN 2739 


CD 


O 

ro 


NATIONAL ADVISORY COMMITTEE 
FOR AERONAUTICS 

TECHNICAL NOTE 2739 



NUMERICAL DETERMINATION OF INDICIAL LIFT AND MOMENT 

FUNCTIONS FOR A TWO-DIMENSIONAL SINKING 

AND PITCHING AIRFOIL AT MACH 

NUMBERS 0.5 AND 0.6 

By Bernard Mazelsky and Joseph A. Drischler 

Langley Aeronautical Laboratory 
Langley Field, Va. 



Tt 


11^ It i * 




AFL 2811 


TECH LIBRARY KAFB, NM 




IG 


MTIOML AI3VISORY COMMITTEE FOR AERONAUTICS 


TECH LIBRARY KAFB, NM 



■ •111 ll.ll Itl.l .lUi III 

ODLSTll 


TECHNICAL NOTE 2739 


NUMERICAL DETERMINATION OF INDICIAL LIFT AND MOMENT 
FUNCTIONS .FOR A TWO-DIMENSIONAL SINKING 
AND PITCHING AIRFOIL AT MACH 
NUMBERS 0.5 AND 0.6 

By Bernard Mazelsky and Joseph A. Drischler 

SUMMARY 


The indicial lift and mcsnent functions are determined approximately 
for sinking and pitching motion at Mach numbers M of 0.5 and 0.6. 

These functions are determined from a knowledge of the existing oscil- 
latory coefficients at the low reduced frequencies and from approximate 
expressions of these coefficients at the high reduced frequencies. 

The beginning portion of the Indicial lift function associated 
with an airfoil penetrating a sharp-edge gust in subsonic flow is 
evaluated by use of an exact method. By use of an approximate method 
for determining the remaining poz*tlon, the complete Indicial gust func- 
tion is determined for M = 0.5, M = 0. 6, and M = 0.7. 

All the indicial lift and moment functions are approximated by an 
exponential series; the coefficients which appear in the exponential 
approximations for each Indicial function are tabulated for M = 0.5, 

M = 0.6, and M = 0. 7. 


INTRODUCTION 


In the study of transient flows, two types of airfoil motions have 
had special significance - a harmonically oscillating airfoil and an 
airfoil experiencing a sudden change in angle of attack. The lift and 
moment functions for an airfoil experiencing a sudden change in angle 
of attack are commonly referred to as indicial lift and moment func- 
tions. In references 1 and 2 the indicial lift and moment functions 
of a two-dimensional airfoil for sinking and pitching motion in a sub- 
sonic compressible flow are presented for M = 0.7. The functions are 
determined from the reciprocal relations between the indicial functions 
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and the oscillatory coefficients. Examination of the oscillatory coef- 
ficients calculated by Dietze, Schade, and P'osslo (see, for example, 
refs. 3^ and 5) indicates that data are not available over a suffi- 

cient range of reduced frequencies for the direct utilization of this 
method to determine the Indlclal lift functions at other subsonic Mach 
numbers . 

The indicial functions for M = 0.8 have been evaluated by Lomax, 
Heaslet; and Sluder, in reference '6, employing an alternate method 
applicable for both subsonic and. supersonic Mach numbers, The alter- 
nate method for subsonic Mach numbers is characterized by an exact, 
readily evaluated expression for the beginning portion of the Indicial 
functions and an approximate expression- for the remaining portion. The 
solution for the remaining portion of the indicial function is lengthy 
and tedious, and, consequently, an alternate approach is used in the 
present paper to calculate this portion of the indicial functions. 

In the present paper indicial functions for M = 0.5 and M = 0.6 
are obtained in approximate f*orra for an extensive range of airfoil 
travel. The beginning portions of the Indicial functions are calcu- 
lated by the method of reference 6. The difficulties of calculating 
the remaining portions of the indicial fimctlons by the method of 
reference 6 are avoided by approximating the missing portions of the 
oscillatory coefficients of references 3 and i)- so that the method of 
references 1 and 2 could be used. The approximation is made in terms 
of the initial parts of the indicial functions and the known parts of 
the oscillatory coefficients. An estimate of the reliability of the 
indicial functions is also presented. 

During the last stage of preparing the paper, it was found that 
Tlmman, Van de Vooren, and Greldanuis ,ha-ve computed the oscillatory 
coefficients in terms of known functions for Mach numbers 0.35, 0.5^ 

0 . 6 , 0 . 7 ^ and 0.8 and for reduced frequencies ranging from 0.1 to 3*0 
(see ref. j) . A discussion of the effects of differences between the 
data of reference 7 and that of the present investigation on the results 
of the calculations is included. 


SYMBOLS 


s distance traveled, half-chords 

0 ) angular frequency 

V forward velocity of airfoil 


c 


chord 
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x,z,^jT[ Cartesian coordinates 

k reduced- frequency parameter, coc/2V 

kQ highest reduced -frequency for whlcii oscillatory 

coefficients are known 


Kko 


^00 


^max 


L(s) 

M(s) 


P 

h 

e 

q 


ki(s) 


mi(s) 



value of reduced frequency where maximim error in 
approximated oscillatory coefficients occurs 

error in oscillatory coefficients at infinite reduced 
frequency 

maximum error in oscillatory coefficients for 
region k > 1 

section lift, positive in downward direction 

section pitching moment about quarter-chord point, 
positive when moment tends to depress trailing 
edge 

dens ity 

amplitude of vertical or sinking displacement, 
positive downward, half -chords 

angle of pitch, positive when trailing edge is 
lower than leading edge 

rate of change of angle of pitch with respect to 
distance traveled in half-chords, positive when 
trailing edge is falling with reference to leading 
edge 

Indlcial lift function for a two-dimensional air- 
foil experiencing a sudden change in vertical 
velocity 

indicial moment function for a two-dimensional air- 
foil experiencing a sudden change in vertical 
velocity where moment is taken about quarter - 
chord point 

indlcial lift function for a two-dimensional air- 
foil experiencing a sudden change in pitching 
velocity about its leading edge 
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M 

Cp 

%(s) 

U ■ 

^ 2 ( 3 ) 

Fc(k) + iGc(k) 

M(k) + iW(k) 

Fcq(k) + iGcq(k) 

Mq(k) + iNq(k) 

f(k) = Fc(k) - F^i 
m(k) = M(k) - M(oo; 
fq(k) = Fcq(k) - I 

mq(k) = Mq(k) - M, 
A,B 

Subscript; 

3c/4 


Mach number 

pressure coefficient . / , ..7 

indicial moment function for a two-dimensional air- 
foil experiencing a sudden change in pitching 
velocity about its leading edge, moment taken 
about quarter-chord point. .. ... 

velocity of sharp-edge .gust 

indicial lift function for a two-dimensional air- _ 
foil due to penetration of a sharp-edge gust 

complex compressible os'cillatory lift derivative 
due to sinking velocity 

complex compressible oscillatory moment derivative, 
due to sinking velocity 

complex compressible oscillatory lift derivative 
due to pitching velocity only 

complex compressible oscillatory moment derivative 
due to pitching velocity only 



( “) 

constants in exponential approximations 
constants defined by equations ( 13 ) and (l4) 

I . 

for airfoil rotating about three-quarter-chord point 
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METHOD OF AMLYSIS 


The indie ial functions are determined in two p^ts. In the region 

0 < s ^ readily detemined indie ial functions given in refer- 

- - 1 + 2M - 

ence 6 are used. For the region s > the relatively more dif- 

X "f" 

ficult method of reference 6 is avoided by an alternate approach which 
utilizes the reciprocal relations between the indie ial and the oscil- 
latory functions given in references 1 and 2. Before these relations 
could be applied at M = 0.5 and M = 0 . 6 , it was necessary to deter- 
mine the unknown oscillatory coefficients at the hi^ reduced frequen- 
cies, k > 1. Approximate expressions for these coefficients are 
obtainable by combining the existing data for the oscillatory coeffi- 
cients at the low reduced frequencies, 0 ^ k 5 1, with the indicial 

functions for the region 0 ^ s ^ ' through the -use of reciprocal 

1 + M 

relations. The reciprocal relations are then used to obtain the 
remaining portion of the indicial functions. 

The following procedure was used to determine the oscillatory coef' 
flcients at the hi^ reduced frequencies: The reciprocal relation 

between the indicial lift function due to sinking velocity and the in- 
phase oscillatory lift coefficient is (see eq. (l4) of ref. l) 


ki(s) = Fc(“) + I r f(k)sin ^ ^ Q) 

K Jq ^ 

or alternately, 

rm.in ks ^ ^ I _ g (3^^., 

JkQ k ^ 2 Jq k 

where f(k) = F^Ck) - Fc(«). 

In equation (lb) values of k 3 _(s) are known from the method of 
reference 6 for 0 < s < value of Ff>(«>) which has been 

shown in reference 1 to correspond to the value of k]_(s) at s = 0 
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is, consequently, also known. The oscillatory coefficient f(k) is 
considered- known for the range of reduced frequencies 0 < k 5 kQ. 
Thus, the ri^t-hand side of equation (ih) is known for values of 
2 M 

0 ^'s ^ In the region kn < k < 00, f(k) 

1 + M 

to be determined. 


is unknown and is 


A solution for. f(k) in the region ■ ko < k < 00 can-be expressed 
by a finite convergent series with imknown coefficients. Definition 
of th^se coefficients is possible by a simultaneous solution of a set 
of algebraic equations. The set -Of eqixations is formulated by sub- 
stituting the finite series for f(k) (together with appropriate 
boundary conditions) in the integral on the left side of equa,tion (lb). 
This Integral is evaluated to obtain an algebraic .expression in terms 
of the unknown series coefficients and the variable s.‘ Each simul- 
taneous equation in the set is then obtained by assigning to the 
variable 3 a value selected in the range wHere ki(s) is known. 

The number of these eqimtions should correspond to' the number. of 
series coefficients to be determined. The set of equations is then 
solved simultaneously for the series coefficients in the approximation 
for f(k). 

The accuracy of the approximate solution for f(k) depends on 
the range of s over which ki(s) is known, the value of kQ, and 
the number of tems in the series. The error~introduced by a finite 
series for f(k) is minimized by the fact that the oscillatory coef- 
ficient for the higher reduced frequencies is associated primarily 
with the Initial portion of the indlclal function which is known. 

Once the approximate expression for f(k) is detemined at the 
high reduced frequencies, k > 1, this expresslpn and the oscillatory 
coefficients at the low reduced frequencies, ' 0 ^ k < 1, can be sub- 
stituted into equation (la) to determine kp(s) for values of 
2 M - ■ 


s > 


1 + M 


The same over -all method may be used for "determining the indlclal 
moment on a sinking airfoil as well as the indlclal lift and moment 
functions' for a pitching airfoil. The reciprocal relation used for . 
approximating the oscillatory moment coefficients at the high reduced 
frequencies due to sinking motion is as fo-llows: 


mQ_(s ) = M(oo) + — 

^ I 


m(k)sin ks 


dk 


(s > 0) (2) 
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where m(k) Is defined 


m(k) = M(k) - M(oo) 


For the case of a pitching airfoil, it was shown in reference 2 that 
the only additional indie lal function that need be considered is the 
one associated with the pitching velocity of the airfoil. The recip- 
rocal relations used in this analysis for determining the oscillatory 
lift and moment coefficients at the hi^ reduced frequencies as well as 
the complete Indiclal functions for the case of pitching velocity are, 
for the lift. 


1 f \ -n. / ^ 2 r"" "kB » / V 

^l„(s) = Fc (cc) + ^ / dk (s > 0) (3) 

'1 Q. ^ ^ 


where 




(i) 


and, for the moment. 


m^(s) = Mq(oo) +1^ £ (k)sln ks ^ 


where 


m^(k) = Mq(k) - Mq(oo) 


( 6 ) 


For convenience, a summary of the reciprocal relations for lift and 
moment due to both sinking and pitching motion is given in table I 
together with the definition of the functions appearing in the 
reciprocal relations. 
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AVAILABLE OSCILLATORY COEEFICIENTS AND THEIR' APPROXIMATION 
AT . HIGH REDUCED FREQUENCIES 


Summary of available- oaclllatory coefficients- at M = 0.^ and 
M = 0.6 .'- The tvo- dimensional compressible osLclUatory coefficients 
for the real and imaginary parts of the lift, and moment for sinking 
and pitching motion are. given in table II -for M = 0.5 and M = 0.6. 
The results given in table II were taken from' two sources (see refs. 3 
and 4) and have been converted to the form given by the equations 
noted in the last column of table I for the JLift and moment. The 
range of reduced frequencies taken from each source is indicated in 
table II. ■ Expressions for. the end points ofrthe oscillatory coeffi- 
cients, that is, for k = 0 and k = », were determined from the 
expressions given in references 1 and 2. 

■ Expressions for indie lal lift and moment functions in region ' 

0 < s < - The expressions for the indiclal lift and moment func- 

~ 1 + M ■ ■ 

tions due to sinking and pitching- motion that are given .in reference 6 

2M 

for the region 0 < s < have been converted to the notation 

" " 1 + M 

used in this paper and are as follows: 


ki(s) = 




- M) 


( 7 ) 


mx(s) = - 


2 




(1 . „) . ^(2 . M) 


( 8 ) 


= m 


- 1 ) 


(9) 


™lq(s) 


. iJ -1 . |£(1 

rtM 1 12 o M 


M) + -(2 

32M^ 


6M + 3m2) 





( 10 ) 
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Approximation of oscillatory coefficients at hi^ reduced fre- 
quencies . - The following series was assumed as an approximate expres- 
sion for f(k) in the range kQ < k ^ 


f(k) = XI anke^^^"^^ (ll) 

n=2 - • 


where aj^ and oc^ axe the unknown coefficients to be determined for 

f(k) at a given .Mach nimiber. Although several forms of series could 
be used to express f(k), this form was fomd to be the most convenient 
for purposes of integration. Note that equation (ll) satisfies the 
condition that 

lim f(k) = 11m |ic(k) - Fc(«>^= 0 (12) 

k^^^ CO k ^ 00 


} 


Two boundary conditions which require the approximate solution 
for f(k) in the region kQ ^ k ^ co to Join the known value at k = kQ 
continuously in magnitude and slope a^e 



(13) 





These two conditions are used to reduce the number of unknowns in 
equation (ll). Substitution of these two conditions into equation (ll) 
leads to the following expressions for the coefficients a]_ and a 2 '. ' 


and 


ai 


B + Axi2 
cxg - aq 


(15) 


9.2 ~ 


gqA + B 
®1 “ °2 


(16) 


I 
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If the expression for f(k) given hy equation (ll) is substituted 
into the left-hand side of equation (ih) and if the values of aq and 
a£ are replaced hy the expressions given by equations (15) and (16), 
then two expressions from which aq and 02 can be evaluated are 
obtained. For the case M = 0.5 and M = 0.6, the value of kQ 
unity, and the two desired expressions are 


B + Aog sq cos sq + aq sin sq 

“2 - °^1 aq2 + sq2 


aqA + B s'-j^ cos s-]_ + a2 sin Sq 
ai - 0-2 .. . 0,22 + sq2 




f(k)sin ksq 

dk 

k 


(iTa) 


and 


B + Aog 82 cos S2 + aq sin 

^2 ~ “'1 ttq ^ + S 2 ^ 



aqA + B S2 cos S2 + CX2 sin S2 
“"I “ °'2 . . ^ 2 ^ 



f(k)sln ks2 


dk 


(1Tb) 


where an and an are numbers greater than zero, and sn and So 

- ■ 2M ^ 

are greater than zero but less than or equal to . 

- 1 + M 

A' simultaneous solution of the nonlinear equations (iTa) and (l7b) 
for the coefficients aq and a2 was obtained by a graphical procedure 

In some cases only an approximate solution for aq and 02 could be 
obtained. In these Instances, values for aq and 02 were selected 
such that the magnitudes of error Incurred in the indiclal function 
were equal at the points Sq and S2» For“the case of M = 0.5; values 

of sq = 0. lOn and S2 = .0.20rt were chosen, and, for M = 0.6, 

Sq = 0.12rt and S 2 = 0.23« were chosen. Both of these values are 

2M 

within. the limits 0 <s < and are approximately equally spaced. 

” “ 1 + M 
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Values for ttie constants a^ and a 2 could then be determined from 
equations ( 15 ) and (I 6 ). In table III values for the constants a^_, 
°'1 j which appear in the exponential approximation 

(eq. ( 11 )) are summarized for the lift and moment due to both sinking 
and pitching motion. 


RESULTS AND DISCIBSION 


Determination of indie ial functions . - With the coefficients in 
equation (11) deteimined for the lift and moment oscillatory coeffi- 
cients due to both sinking and pitching motion, the indicial functions 


(s) were evaluated over the range 


k]^(s), m 2 ^(s), k]_ (s), and m-| 

2M ^ 

s > T J 7 ’ for M = 0.5 and M = 0.6 by the methods of references 1 

1 + M ' 

and 2. These results, and corresponding functions calculated by the 

< 2M 


method of reference 6 for 0 < 


s i 


are shown, respectively, in 


= 1 + M 

figures 1, 2, 3j and 4 together with the results for M = 0,7; obtained 
from references 1 and 2. Also included is the function for incompressi- 
ble flow obtained by Wagner. The indicial lift and moment functions due 
to pitching velocity for the airfoil rotating about the three-quarter- 
chord point are shown in figures 5 and 6. For incompressible flow, . 
the lift and moment are impulsive at s = 0 and constant thereafter; 
however, for compressible flow, the lift and moment are finite and 
time-dependent and approach a constant value very rapidly. Comparison 
of the Indicial lift functions due to both sinking and pitching 
velocities for subsonic Mach numbers, figures 1 and 3; indicates that 
as the subsonic Mach number increases the lift at s = 0 becomes less 
and, for s equal to more than approximately two half-chords, the 
growth of lift to the steady state becomes less rapid. For the case 
of the moments about the quarter- chord point due to si nking and 
pitching motion, it is similarly seen from figures 2 and 4 that while 
the moments are impulsive at zero for incompressible flow the moments 
are finite and their values decrease with increasing Mach number for 
compressible flow. Although the moments for s > 0 are constant for 
incompressible flow, a time -dependent moment which rapidly approaches 
a constant exists for compressible flow where- the rapidity of approach 
decreases with Increasing Mach number. 


Considerations for indicial lift due to penetration of a sharp-edge 
gust .- In references 1 and 6 it was shown that the indicial lift func- 

2M 

tlon k]_(s) for the region s > is associated primarily with 

1 - M 

the lift due to circulation only. As a consequence, an approximation 
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for the indicial lift function k 2 (s) due ; to penetration of a sharp- 
edge gust in this region can he determined on the basis "of the relation 
ship that exists between the two indicial lift functions _ ki(s) and 
kgCs) for the circulatory portion of the lift for incompressible flow. 

In reference 1 the function k 2 (s) is determined for the region 

2M ' ' 

3 > ^ on the basis of this relationship for M = 0.7* However, 

no information concerning the beginning portion of this function is 
available at present, 

A method for determining an explicit expression for this function 

in the region 0 5 s ^ ^ is possible from an extension of the 

analysis given in reference 6 where the boundary conditions for the 
perturbation velocity due to the penetration of a gust can be sub- ’’J 
stituted for the perturbation velocity due to sinking motion. Details 
of the evaluation of the function k 2 (s) for this region are given in 

appendix A. The expression obtained -for this function is relatively 
simple as can be seen from the following equation: 

( 18 ) 


where k 2 (s) is defined such that the lift on the airfoil is 


L(s) = -npcVUkgCs) 


(19) 


The gust function k 2 (s) 


for the region 


_^<s <-^ 
1 + M= =1-M 


can 


be determined from an extension of the results given in reference 6. 
An explicit expression cannot be easily obtained because of the com- 
plexity of the downwash regions that must be considered. Since the 
penetration function is known everywhere except in the region 


2M 


2M 


< s ^ a reasonable estimate of the curve may be obtained 

- - - 1 - m' •'2M 

and 


1 + M 

by joining the ends of the two solutions , at the points 
2M 

by a smooth curve. 

1 - M 


1 + M 


In figure 7^ the Indicial lift functions due to penetration of a 
sharp-edge gust are shown for M = 0.5, M = 0.6, and M = 0.7 ^or an 
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extensive range of airfoil travel. Also shown is the solution for M = 0 
obtained by Kiissner. Inspection of figure 7 indicates that not only is 
the growth of lift to the steady state less rapid as the subsonic Mach 
number is increased but also the magnitude of lift in the beginning 
portion is less. 

Approximation of indicial functions by analytic expressions .- Since 
the exponential function has a simple operational equivalent and a series 
of such functions has been found convenient for approximating k]_(s) at 
M = 0 and M = 0.7^ a finite sum of exponential terms Was chosen to 
approximate the Indicial lift and moment functions for sinking and 
pitching motion at M = 0.5 and M = 0.6 as follows: 

bo + b]_e“^l® + b 26 + b^e ^3® (20) 


In table IV values for these constants are given for the lift and 
moment due to sinking motion as well as for the lift and moment due to 
pitching motion about the three -quarter-chord point. The constants 
are noted for M = 0.5 and M = 0.6 together with those for M = 0.7 
from reference 2. In reference 8 the constants are given for incom- 
pressible flow. Plots of these exponential approximations are shown 
in figures 1, 2, 5^ and 6. Comparison of the approximations with the 
actual curves in these figures indicates good agreement. 

The exponential approximation given by equation (20) was also 
found to be convenient for fitting the indicial lift functions k 2 (s) 
due to penetration of a sharp-edge gust shown in figure J. In table IV 
the constants required for fitting the’ k 2 (s) function at M = 0.5; 

M = 0.6, and M = 0.7 for use in equation (20) are given. In refer- 
ence 8 the constants for this indicial f met ion are given for incom- 
pressible flow. Comparison of these approximations with the derived 
curves (see fig. 7) Indicates excellent agreement. 

Approximation of oscillatory coefficients by analytic expressions . - 
The corresponding approximate expressions for the harmonically oscil- 
lating airfoil can be fomd from equation (20) and an alternate form of 
the reciprocal equations. For the case of the lift due to a i nking 
motion, the reciprocal equation is 


Fc(k) + iCc(k) = ik 



e 


-iks 


ds 


( 21 ) 
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Substitution of the' approximation for k]_(s) given by equation (20) 
leads to the following expressions for 5'c(k) acd Gc(k) 


F^(k) = bo + 


b-i k^ bpk^"' bok^ 

+ — r — 


+ k^ 02 ^ + ^ 33^ + 


> ( 22 ) 


b]_|3j_k 1 > 2 ^ 2 ^ -^ 3 ^ 3 ^ 

Gc(k) = + — r r- + 


3^^ + k^ 02 ^ + k^ 100^ + k^ 


In a similar fashion, expressions for the moment functions M(k) and 
W(k) can be obtained in the form given by equation (22) as well as 
for the lift fimctions due to pitching velocity about the three- 
quarter-chord point moment 

functions (Mq) 2 Qyi^(k) and ' (Wq) 2 c/ij.(k), The transfer 'of these coef- 
ficients from the leading edge to the three-quarter-chord point was 
made with the aid of- equations ( 27 ) of reference 2. 


Estimate of reliability of indicrial 'functions . - Although the 

2M 


Indicial functions are known exactly in the region 0 < s < 


= 1 + 


the 


reliability of these functions for values of s > is question- 

1 + M . - 

able since the oscillatory coefficients at the. high reduced frequencies, 
which were used in determining the indicial function in this region, 
have been determined only approximately. An estimate of the reliability 
of the Indicial functions in this region due to the errors in the oscil- 
latory coefficients at the high reduced frequencies could be made on 
the basis of the differences which exist between the approximate solu- 
tion (based on the oscillatory coefficients) and the exact solution 

(based on ref. 6) of the indicial functions in the region 0 < s < — . 

1 + M 

Since the oscillatory coefficients at the high reduced frequencies 
principally affect the beginning portion of the function k]_(s), the 
differences between the approximate and exact ^solutions in this region 
should give a reasonable and conservative estimate of the errors 
Incurred in the indicial functions for larger values of s. 
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In appendix B, a chart (see fig. 8) is obtained for indicating the 
variation of errors in the indicial functions with airfoil travel s 
due to the errors in the in-phase oscillatory coefficients at the high 
reduced frequencies. Although figure 8 indicates the errors in only 
the indicial lift function k]_(s), the results in this figure are also 
applicable for the other indicial functions. Inspection of figure 8 
for kQ = 1 indicates that the maximum error in the indicial function 
occurs within q. quarter of a chord regardless of the. position of the 
maximum error occurring in the in-phase oscillatory coefficients (as 
Indicated by several values of k). Since s = 0. 5 lies within the 
PM 

region 0 5 s ^ for M = 0. 5 and M = 0.6, it may be expected 

1 + M ■ 

that the greatest errors in the indicial functions occur in this region 
and that errors beyond this point approach zero very rapidly. The 
largest differences between the Indicial functions determined from the 
oscillatory coefficients and by the method given, in reference 6 in the 
2M 

region 0 < s < — were found to be less- than 1. 5 percent for 

1 + M 

M = 0.5 and M = 0.6. Thiis, it appears from figure 8 that the largest 

2M 

errors in the indicial functions beyond the point s = are much 

1 + M 

less than 1. 5 percent and appear to approach zero very rapidly for 
increasing values of airfoil travel s. 

Effects of discrepancies between data of reference 7 and those of 
present investigation on indicial functions . - Just prior to completion 
of this paper oscillatory lift and moment coefficients which cover a 
larger range of reduced frequencies than that considered by Dietze and 
Schade were published (see ref. 7)« 1°. figures 9 and 10 a comparison 

is shown of the in-phase coefficients used in the present analysis 
with those obtained from reference 7- For the range of reduced fre- 
quencies 0 < k ^ 2.5j the data agree within approximately 3 percent. 
However, for k > 2.5 the differences tend to become larger than 
3 percent in some cases. 

Consideration of the over-all effect of the discrepancies shown 
in figures 9 and 10 indicates that these differences would result in 
a maximimi error of 3 percent' in the indicial functions. This amount 
of error does not appear serious, and in view of the few data avail- 
able from reference J, recomputation of the Indicial function does 
not appear worthwhile. It might be noted that the data of Dietze and 
Schade agree better mutually than with the data of reference 7- The 
reason for these discrepancies is not known; consequently, the selection 
of the most accurate oscillatory coefficients is still questionable. 
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CONCLUDING REMARKS 


From approximate expressions for the in-phase oscillatory coeffi- 
cients determined for the high reduced' frequencies and with a knowledge" 
of the existing data at the low reduced frequencies, the indicial lift 
and moment functions have been determined for sinking and pitching 
motion at Mach numbers M of 0.5 and 0.6. ’ 

The indicial lift function due to penetration of a sharp-edge gust 
is determined for M = 0.5, M = 0.6, and M = O.J- for an extensive 
range of airfoil travel. Inspection of these functions indicates that 
not only is the growth of. lift to the steady state less rapid as *the 
subsonic Mach number is increased but also the magnitude of lift in 
the beginning portion is less. ‘ 


Langley Aeronautical.- Laboratory, 

National Advisory Committee for Aeronautics, 
Langley Field, Va. , April 15, 1952. 
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APPENDIX A 


DETERMIMTION OF INDICIAL LIFT FUNCTION k2(s) DUE TO PENETRATION 

OF A SHARP-EDGE GUST FOR 0 < s < — 

- 1 + M 


In reference 6 a mettiod is outlined for determining the suhsonic 
indicial lift functions for a sinking airfoil. However, if the boundary- 
conditions for the pertiirbation velocity due to the penetration .of a 
gust are substituted for the perturbation velocity due to sinking motion, 
the function k 2 (s) can be obtained. Because of the complexity of the 
problem arising after the Mach forecone intersects the trailing edge, 

2M 

the analysis is made only for the region 0 = s = The boundary 

condition for the gust velocity U differs from the boundary condition 
of the sinking airfoil in that there is no perturbation velocity on the 
part of the airfoil left unshaded in the following sketch: 
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If as in. reference 6 the axes in the preceding sketch are transformed 
by the relat,ions , v / 



then the geometry can be .represented by 



The expression for the pressure coefflclerjt ' Cp can be written in 
terms of the coordinate system as follows, where the subscripts 

indicate the region for- which the particular equation applies: 

Cp^ = 0 (Al) 



hij J ^ /I . /t, - Tl^ 


(A2) 


Here it might be noted that equation (A2) differs from the corresponding 
equation in reference 6 only by the lower limit of the last integral 
since this integration extends over only that part of the wing which . 
has penetrated the gust. 
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Evaluation of equations (Al) and (A2) gives 


Cpi = 0 


(A3) 


Jm - ^ 


^P2 = 


rt(l + M)l/^ + X 


(A4) 


and the expression for the Indlclal function k 2 (s) "becomes 


’sc 


/2M 


■(C - 


sc 


2ick2(s) = i 


C-D dx + — / C-n, dx 

-sc/2 . 2 sc/2M 


fo<s<-^) 

1 + v 


or 


(A5) 


k2(s) = 

rt/M 


^ ^ 8 < 


2M 


1 + M. 


(A6) 
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APPENDIX B ■ 

INVESTIGATION OF ERRORS IF INDIO lAL FUNCTION DUE TO ERRORS IN IN-PHASE 
OSCILLATORY COEFFICIENTS AT HIGH REDUCED FREQUENCIES 

Consider the following reciprocal relation for the lift due to 
sinking velocity; 

k, (s) = 2 r Fcfr)sln to (Bl) 

^ 0 ^ 

It is assumed that F(,(k) is known for values of 0 k. ko; however, 
for values of k > kQ an error in the F^(k) function is assumed of 
the following fonn: 
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If equation (B2) is substituted into equation (Bl), an expression for 
the error in the ki(s) function may be obtained as follows: 



= 


^ - sin kQS - Si(koB) + kQS Ci(kos) 


(B3) 


where &i(x) and Ci(x) are tabulated in Jahnke and Emde (ref. 9)» 

If the value of is known and if the maximum error is 

assumed to occur at some arbitrary point nkQ (where K > 1), two 
error functions of the type given by equation (B2), where one is dis- 
placed from the point kQ by an amount KkQ, may be superposed such 

that the value of = 0 as Illustrated in the following 

sketch: , ^ 

Fc(k) 



The value of the maximum error which occurs at k = kIsq can be 
expressed in terms of €00 as follows: 


^max 
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The corresponding expression for the error in the function k2_(s) due 
to the total error in the function shown in the preceding sketch 

is 


2s, 


max 


kx(s) 


error 


■(^l^s Cl(kQs) - SKko Ci(Kk()s) + Si(KkQs)'- 


Si(kos) + sin(Kkos)-- sin(kQS^ (B 5 ) 

A plot of equation (B 5 ) as a function of Icqs is shown in figure 8 
for several values of the parameter k. Thus^ for a maximum error %ax 
in the Fc(k) function at any point ,ko < «ko < an estimate of the 
corresponding error in the function ki(s) can be obtained from 
figure 8 . 


r 


( 
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>PATffg 1 

SUM4AAY OF RBC1FH0CAL RBLAl^ONa TOR LIFT AHD MCMERf StlE TO dlnOJO 
ASD PtZCfflBa MCTIOH UT SUffiOinC FLOU. 


DAflnltlon of IsdlclAll 

I of oeclllatory ooofflciaat' 




Lift 


t (.) . 1 ^( 0 ) * i /’* SgW°°* 4t 

* Jq It 

• « 2 

L(.) - »i(b) 

L(i) - -apc72a^^(-lhJt)[rpClt) + IflcCki] 


MCQMOt^ 



.i(.) .2 rmiMiaJsLii. 

‘>'0 k 

. K( 0 ) t 2 r »(“)«”■ 4 k 
. *>^o k 



(g) - epc®v^ 5 ^*i(b)| **(■) ■ >tPc8T8e^(-lhk)5(k) + lK(kJ] 




Lift 





ki{B) - rc(0) t ^ 

■ 

1 

1 

• 

L{g) - -icpcT2ee^[Fe(]c) + lOcU]] 


kcBgi^f_ 
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.i('.l at 
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«l(. ) . *(0) t 2 r "(IQ«=°« 2(«) 4k 
rt ^0 

M(g) - H3C^f2ftn2_(a) 

M(«) • «pcar2o«^* j5(k) + i*(k3 


Fca(k)gln In 


w / X -a /’•^caWg 


5 • , , 2 /^OcoWco* !“ I 

f 


L(e) ■ .«pcV*22lti^(») !,(•) . -iiJKV221Wd‘'“^5^(lt) + 



V‘>. 


H(>) - HJc&I?2yil^(«) M(l) . «pc2v22iU4“’[^(k) + l»j{lc^ 


^Hocgat t«]nQ a'bout quartar-cbord point 




i'H... n'’nt 7 ^-aiiT 























TABLE II 


oaciujuraRY coeticiehis for oeb lift asd mombbt dub to simaiB Am 

PUCHIHG MMICIf AT H = 0.5AroM = 0.6 


Q 


Reference 

■ 

Fo(k) 

for - 

- 3 c <k) 

for - 


for - 

-0cq(4) 

for - 

HI 


-ir(k) 

for - 


for - 

-H^(lc) for - 



















■ 

B 




B 

M = 0.6 

B 

M » 0.6 


H 0 0.6 


0.6 



BB 

H > 0.6 

I and 2 

0 

IQQ 

1.2300 

■ 

H 

0.8660 

0.9375 

■ 

0 

0 

0 

0 

0 

0.07217 

0.07812 

m 

■ 


.02 

1.0962 

1.1747 

.1042 

.1307 

.8323 

.ffr 98 

.09675 

..1177 

.000125 

.00025 

.00375 

.004500 

.07327 

.08397 

.002625 

,002250 


,04 

1.0385 

1,1010 

.1517 

.1885 

.7777 

.8239 

.1191 

.1496 

.00050 

.00100 

,007375 

.009125 

.07292 

.07890 

.004825 

.0058^5 


,o6 

.9858 

1. 0336 

.1766 

.2176 

.7362 

.7768 

.1413 

.1756 

.0008333 

.001667 

.01092 

.01317 

.07200 

.08037 

.006619 

.007967 


.08 

.9394 

. 9 ^ 

.1885 

.2363 

.7032 

.7728 

.1547 

.1614 

.001312 

.002500 

.01425 

.01712 

.07347 

.08075 

,009009 

.01054 


.10 

.8980 

.9313 

.1919 

.2321 

.6732 

.6938 

.1657 

.1972 

.001700 

.003300 

.01760 

.02100 

.07369 

.08083 

.01070 

.01254 

3 

.20 

.7648 

.7830 

.1535 

.1877 

.5693 

.5782 

.1490 

.1729 

.004075 

.007500 

.03393 

.03875 

.07552 

.08451 

,02020 

.02260 


.30 

.6938 

.709a 

.08966 

.1205 • 

.5155 

• 5201 . 

,U76 

.1463 

.006483 

.01187 

.04907 

.05600 

.07^6 

.08809 

.02912 

.03173 


.40 

.6607 

.6759 

.02475 

.0559 ' 

.4864 

.4900 

.08310 

.1147 

.009250 

.01693 

.06475 

.07312 

.06014 

.09254 

,03780 

.04033 


.50 

.6447 

.6645 

-.03750 

-.002200 

,4708 

.4756 

.05249 

,08808 

.01260 

.02324 

.0B050 

.09020 

.08310 

.09768 

.04652 

.04850 


.60 

.6415 

.667a 

-.09550 

-.05383 

.4644 

.4695 

.02498 

.06646 

.01675 

.03127 

.09667 

.1067 

,08641 

.1042 

.05509 

.05627 


.70 

. 61^79 

.6804 

-.1502 

-.^36 

.4644 

.'4690 

.00041 

.04881 

.02200 

.04159 

• .1130 

.1244 

.09047 

,1105 

,06357 

.06346 

1 , 

.80 

.6619 

.6992 


-.1377 

.4665 

.4707 


.03683 

.02737 

.05116 

.1293 

.1409. 

,09385 

,1165 

.07192 

.06975 


1.00 

.7105 

.7636 


-.1890 

.4861 

,4839 


.02883 

.04405 

.08590 

,1643 

.1727 

.1957 

.1364 

.06864 

,07806 

1 and 2 

ee 

1.2732 

1.0630 

■ 

■ 

.6366 

.5305 

■ 

■ 

.3183 

.2653 

■ 

■ 

.2653 

.2210 

HI 

m 




ro 
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TABLE III 


ro 


COEFTICIEHTS IN EXPONElfriAL APPROXIMATION FOR 

OSCILLATORY LIFT AND MCMEHT FUNCTIONS AT HIGH REDUCED 
FREQUENCIES FOR M = 0.5 AND M = 0.6 


Oscillatory 

function 

aq for - 

32 

for - 

ai 

for - 

02 

for- - 

M = 0.5 

M = 0.6 

M = 0.5 

M = 0.6 

M = 0.5 

M = 0.6 

M = 0.5 

M = 0.6 

f(k) 

-0.003203 

-0. 1566 

0. 5595 

0.liK)75 

0.200 

1.000 

1.67 

3.150 

m(k) 

-.07520 

-.02003 

.3i^95 

*199^ 

.385 

.001 

1.256 

2.056 

fq(k) 

-.06812 

0 

.082W 

.01^665 

.584 

00 

3.155 

1.100 

mq(k) 

. 11575 

. 08110 

.(A375 

. 003572 

1.175 

1.400 

2,560 

26.000 
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TABLE IV 


COEFPICIEHTS IN EXPEESSION Lq + bj^e + b 2 e"^^^ + FOR APPROXIMATING 

THE HTDICIAL LIFT AND MCMENT FUNCTIONS AT M = O. 5 , 0.6, AND O .7 


Indie iai 
function 

Mach 

number 

^0 


^2 

’'3 

Pi 

P2 

^3 


0.5 

1.155 

-0.406 

- 0.249 

0.773 

0. 0754 

0.372 

1.890 

ki(B) 

.6 

1.250 

-.452 

-.630 

.893 

.0646 

.481 

.958 

.7 

I.ij-OO 

-.5096 

-.567 

-.5866 

.0536 

.357 

.902 


.5 

0 

.0557 

- 1.000 

.6263 

2.555 

3.308 

6.09 


.6 

0 

-.100 

-1.502 

1.336 

1.035 

4 .o 4 o 

5. 022 

.7 

0 

-.2425 

.084 

-.069 

.974 

. 668 

.438 


.5 

0 

0 

- 2.68 

2.362 


4.08 

4.90 


.6 

.7 

0 

0 

0 

-.083 

0 

-.293 

-.2653 

.149 

.800 

1.565 

1.345 

2.44 


.5 

-.0721 

-.248 

.522 

-.2879 

1.562 

2.348 

6.605 

(s) 

.6 

-.0781. 

-.077 

.380 

-.2469 

.551 

2.117 

4.138 

V V3c/U- 

.7 

-.0875 

-.00998 

.1079 

-.02920 

.1865 

l.l 4 l 

4.04 


.5 

1.155 

-.450 

-.470 

-.235 

.0716 

.374 

2.165 

kgCs) 

.6 

1.250 

■ -.410 

-.538 

-.302 

.0545 

.257 

1.461 

.7 

1 . 1^00 

-.563 

-.645 

-.192 

.0542 

.3125 

1.474 
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Figure 1.- Comparlaon of indiclal lift functions due to a sudden change 
in sinking velocity or pitching displacement for several siibsonlc 
Mach numbers. 
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Figure 2.- Compariaon of Indiclal moment functions due to a sudden 
change in sinking velocity or pitching displacement for several 
suhsonic Mach numbers. (Moment taken about quarter- chord point.) 
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Figure 3-- Comparison of indie lal 1: 
in pitching Telocity for an alrfi 
for several subsonic Mach numberi 
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Figure 4 .- Comparison of indicial moment fmctions due to a sudden 
change in pitching velocity for an airfoil rotating about its 
leading edge for several subsonic Jfech numbers. (Moment taken 
about quarter-chord point. ) 
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Figure 6.- Comparison of Indiclal moment functions due to a sudden 
n hftn In nr 1’hv on a -T T -rv\4*o-n*«r^ ^4-,, 

- — 1_J JC" —a ' — J. w uq^ u -1-J_i.gj OrU^UU X UD 

I three-quarter-chord, point for several subsonic Mach numbers, 
(Moment taken about quarter-chord point.) 
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Figure 7« - Comparison, of the indicial lift functions due to penetration 
of a sharp-edge gust at several subsonic Mach numbers. 
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Figure 6. ^ Effect of errors in f (K) at high, reduced, frequencies on 

indicial lift function kiCa). 





(b) Pitching velocity. 


Figure 9,- Comparison of In-phase oscillatory" lift coefficients used in 
present analysis with coefficients determined from reference 7, 












(b) Pitching velocity. 

Figure 10.- Comparison of in-phase oscillatory moment coefficients used 
in present analysis with coefficients determined from reference 7 . 


NACA-Langley - 7-24-52 - 1000 





